Study of Possible Proton- Antiproton Hexaquark State in Lattice QCD 
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We have used standard techniques of lattice quantum chrome-dynamics to look for evidence of 



the spin-zero six quark flavour singlet state (J 



0" 



observed by BES Collaboration, and to 



determine the splitting between the mass of the possible proton-antiproton and the mass of two 
protons, its threshold. Using the various interpolating fields in the quenched approximation we find 
indications that for sufficiently light quarks proton-antiproton is slightly above the 2m p threshold. 
The lattice resonance signal of binding observed near the physical and continuum regimes do not 
support the existence of proton-antiproton state as a spin-zero hexaquark state in quenched lattice 
QCD. 
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Recently, the BES Collaboration in Beijing observed 
a near threshold enhancement in the proton-antiproton 
(jrp) mass spectrum from J/ip — > ^ypp radiative decay 
[lj. Fitting the enhancement with an S- wave Breit- 
Wigner resonance function, results in peak mass at M = 
1859+ 3 (stat)+ 5 (sys) with a total width T < 30 MeV/c 2 . 
With a P-wave fit, the peak mass is very close to the 
threshold, M = 1876.4 ± 0.9 MeV and the total width is 
very narrow, T = 4.6 ± 1.8 MeV. This discovery is subse- 
quently confirmed by the Belle Collaboration in different 
reactions of the decays B + — > K + pp Q and B° — ► D°pp 
[H, showing enhancements in the pp invariant mass dis- 
tribution near 2m p . Such a mass and width does not 
match that of any known particle 

Theoretical existence of possible proton-antiproton 
bound state has long been speculated in quark model and 
conventional nucleon potentials @, @, However, it 

was only in a very recent study, made by Yan et al us- 
ing the Skyrme model [9j, that predicted mass and the 
width close with the experiment. Experimentally, the 
quantum numbers of pp(1857) are not well determined 
yet. The photon polar angle distribution was found con- 
sistent with l + cos 2 # 7 suggesting the angular momentum 
is very likely to be J = 0. Making full use of general sym- 
metry requirement and available experimental informa- 
tion the corresponding spin and parity are J PC = 
[Toj ] . In this letter we report first quenched lattice QCD 
calculation capable of studying the pp(1859). 

In contrast with conventional baryons and mesons, it 
is difficult to deal with q" 1 ^ (m + n> 3) states in lattice 
QCD. For example, a q 3 q 3 state can be decomposed into 
a couple of colour singlets states even in the absence of 
unquenched effect. The two-point function, in general, 
couples not only to the single hadron but also to the 
two-hadron states. We seek operators that have a lit- 
tle overlap with the hadronic two-body states in order to 
identify the signal of our hexaquark state in lattice QCD. 
For this purpose, we first consider a non two-body-type 
field for I(J PC ) = 0(0 and construct our interpolat- 
ing fields based on the diquark description. 



In Jaffe model 11 



each pair of [ud] form a diquark 
which transforms like a spin singlet (l s ), colour anti- 
triplet (3 C ), and the flavour anti-triplet (3/). Therefore, 
for a diquark operator, one has [l2| 



'y 



(x) = ^e i:jk e abc qj tb (x)Crq kyC (x), (1) 



where e a bc is completely antisymmetric tensor, and 
(a, 6, c) and (i,j,k) denote the colour and flavour in- 
dices, respectively. The superscript T denotes the trans- 
pose of the Dirac spinor and C is the charge conju- 
gation matrix. The Dirac structure of the operator is 
represented by the matrices T, satisfying T a fj = — Tp a 
{a and (3 are Dirac indices) such that the diquark op- 
erator transforms like a scalar or pseudoscalar. The 
colour and flavour antisymmetries restrict the possible 
T's within 1, 75 and 757^,. Then the hexaquark hadron 
pp([u(i][ud][ww]) emerges as a member with S = and 
I = in (Q 3 ® Q 3 ) ® Q 6 = ([15 cs ] g> [15 cs ]) <g> [21 cs ] in 
SU(6) colour-spin representation and a flavour singlet in 
(3/(g)3/)®6/. With this picture the interpolating oper- 
ator for ^ is obtained as 



Xd q {x) 



^abc 



1 . a 



(x)Q 



r' 



b 'x)Q'r,f(x) 



(2) 



where Q = u T C 1 7sT u. This identification looks fa- 
miliar if we represent one of the quarks by charge conju- 
gate field: q a q b — > q~c a qb, where q c = -iq T o- 2 -f 5 . Then 
the classification of diquark bilinears is analogous to that 
of qq bilinears. We choose F = 1 and T = T =75. 
There are many more possibilities of constructing the 
operator even in the I — channel. In principle testing 
various other interpolating operators for the best overlap 
with *~ state should provide information on the wave 
function of the particle. Our second operator, which we 
will refer to as the iViV-interpolating field, is generalised 
to obtain an isospin 1 = colour-singlet iViV-type hexa- 
quark interpolating field: 



Xnn 



1 

71 



[O n O n - o R o N \ , 



(3) 



2 



where 



ansatz 



O n {x) 



and 



c -abc [(Ua(x)C'y 5 d b (x))u c (x) - (u <-> d)u c (a:)] 

(4) 



■abc [{u a {x)C"/ 5 d[ (x))u c (x) - (u <-> d)u c (a;)] 

(5) 

are the standard interpolating fields for the upper two 
Dirac components of the nucleon and antinucleon op- 
erators, respectively. These definitions take account of 
the symmetries for isospin and nonrelativistic spin rep- 
resentations. In absence of quark annihilation diagrams, 
correlation function is expressed in terms of basis deter- 
mined by direct (36 members) and cross (324 members) 
contributions of quark contractions. From the dispersion 
of the two-point function, the contribution from direct 
and the cross amplitudes can be easily calculated. We 
have taken single quark-antiquark exchange and a dou- 
ble quark-antiquark exchange to be equivalent because 
the N and N operators are summed over all spatial sites 
for s-wave state. A more accurate procedure will be to 
measure the 2x2 correlation matrix of two different in- 
terpolating operators in Eqs. ((2]) and (J3|) and extract 
the effective mass from the eigenvalues using variational 
techniques [HI]. In this study, however, we do not intend 
to pursue this issue further. 

To examine the possible proton-antiproton ground 
state in lattice QCD, we generate quenched configu- 
rations on a 20 3 x 60 lattice (with periodic bound- 
ary conditions in all directions) with tadpole-improved 
anisotropic gluon action 14| in the coupling range of 
3.0 — 4.0. Gauge configurations are generated by a 5- 
hit pseudo heat bath update supplemented by four over- 
relaxation steps. These configurations are then fixed to 
the Coulomb gauge at every 500 sweeps. After discarding 
the initial sweeps, a total of 500 configurations for each 
coupling are accumulated for measurements. Using the 
tadpole-improved clover quark action on the anisotropic 
lattice 115J over the lattice spacing range of 0.26 — 0.39 
fm, we compute the light-quark propagators at six values 
of the hopping parameter n t which cover the quark mass 
region of m u < m q < 2m u . 

To obtain a better overlap with the ground state we 
used iterative smearing of gauge links and the applica- 
tion of fuzzing technique for the fermion fields [161 ]. The 
application of fuzzing for two of the six quarks inside 
the hexaquark state flattens the curvature of the effec- 
tive mass. The largest plateau in the region with small 
errors is obtained with fuzzed u- and d- quarks. We used 
this variant to calculate our correlation functions. 

From the correlation functions we extract the mass (en- 
ergy) by standard y 2 fitting with multi-hyperbolic cosine 



C(t) = y^v4 t cosh(tmj). 



However, to ensure the validity of our results, we com- 
pared them with those obtained as the solution to the 
equation 

C(t + 1) _ cosh[(i + 1 - N t /2)a t m ef f} 



C(t) 



cosh[(i - N t /2)a t m e ff] 



rmn,tmax] for the final analysis is 



for a given C(t + 1)/C(i) at fixed t 

The fitting range [t 
determined by fixing t max and finding a range of t min 
where the ground state mass is stable against t m i n . We 
choose one "best fit" which is insensitive to the fit range, 
has high confidence level and reasonable statistical errors. 
Typical example of the effective mass plot is shown in Fig. 
[TJ The effective mass decreases monotonically, which im- 
plies that excited spectral contributions are greatly re- 
duced. For the XAW~tyP e field, an impressive plateau 
with reasonable statistical errors is observed in the in- 
terval 35 < t < 50, where we expect to achieve single 
state dominance. For the Xdg~tyP e field, the signal is 
noisy at earlier times, and hence we fit in the interval 
45 < t < 55. The best-fit curve to the M N fj data has 
X 2 /N DF = 1.01 and for M dq , we find x 2 /N DF = 0.89. 
Statistical errors of masses are estimated by a single elim- 
ination jackknife method. We kept statistical errors un- 
der control by ensuring that analyzed configuration are 
uncorrelated, which is made possible by separating them 
by as many as 500 sweeps. The statistical uncertainties 
on our hadron masses are typically on the few percent 
level. The masses of the non-strange mesons tt, p as well 
as the nucleon were computed for scale setting and ana- 
lyzing the stability of hexa-quark state, respectively. 




FIG. 1: Effective mass of the J PC — h colour-singlet state 
as a function of t. The data correspond to interpolators \dg 
(open circles) and Xnn (open triangles) at j3 — 3.5 and n t = 
0.2380. 
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The chiral extrapolation to the physical limit is the 
next important issue. From the view point of chiral per- 
turbation theory, data points with smallest should be 
used to capture the chiral log behaviour. Lcinweber et al 
171 ] demonstrated that the chiral extrapolation method 
based upon finite-range regulator leads to extremely ac- 
curate value for the mass of the physical nucleon with 
systematic errors of less than one percent . Unfortunately 
there does not seem to exist such a chiral extrapolation 
technique for multi-quark hadrons beyond a naive one. 
Since quenched spectroscopy is quite reliable for mass 
ratio of stable particles, it is physically even more mo- 
tivated to extrapolate mass ratio instead of mass. This 
allows for the cancellation of systematic errors since the 
hadron states are generated from the same gauge field 
configurations and hence systematic errors are strongly 
correlated. Fig. [2] collects and displays the resulting par- 
ticle mass ratios, M N ^ dq /M p , extrapolated to the phys- 
ical quark mass value using linear and quadratic fits in 



m/i — a + bm n , rtih = a + hm v + cra^. (6) 

For the nucleon we have employed the procedure adopted 
171 ]. The difference between these two extrap- 



in Ref. 

olations gives some information about systematic uncer- 
tainties in the extrapolated quantities. The data for the 
mass ratios M N fq dq /M p behave almost linearly in 
and both the linear and quadratic fits, in Eq. ©, es- 
sentially gave the identical results. The mass ratios ex- 
tracted from the fields \dq and Xnn are very close and 
agree within the errors. It can be seen that mass ratios 
show a weaker dependence on the quark mass. We expect 
the contributions from the uncertainties due to chiral log- 
arithms in the physical limit significantly less dominant 
at our present statistics. 

The question whether the lowest-lying J PC = h 
state extracted is scattering state or a resonance is bet- 
ter resolved by analysing the ratio of the mass differences 
AM = M N fq^ dq — 2Mp (between the candidate low-lying 
hexa-quark state and the free two-particle state) and the 
nucleon mass. The fitted data, illustrated in Fig. [3l 
show that both the field operators give the consistent es- 
timates of the ratio and the masses derived from both the 
operators are consistently higher than the lowest-mass 
two-particle state. The ratios are clearly positive and in- 
crease in magnitude as we approach the physical regime. 
At zero quark mass the mass difference for both the field 
operators is ~ 55(2) MeV with \ 2 /Ndf — 1-05. In other 
words, binding becomes weaker near the physical regime 
with a general trend of positive binding as the zero quark 
mass limit is approached. The chiral uncertainties in the 
physical limit are less than 2%. 

A similar behaviour observed for the mass differences 
between the hexaquark and the two-particle states at 
smaller lattice spacing a s = 0.26 fm is illustrated in Fig. 
31 The ratio obtained from the Xdq operator for the five 
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FIG. 2: Chiral extrapolation of hadron mass ratios, 
M N ff dq /Mp, extracted with the diquark-type (open trian- 
gles) and NN (open circles) interpolating fields at a s — 0.32 
fm. The solid triangles correspond to the data from proton 
interpolating field. The dashed lines are chiral extrapolation 
(linear in m^) to the mass ratios using Eq. (|6]). 
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FIG. 3: Chiral extrapolation of the ratios, AM/M P = 
{M N ff dq — 2Mp)/Mp, extracted from diquark-type (open cir- 
cles) and NN (open triangles) fields at a s = 0.29 fm. The 
solid and the dashed curves are chiral extrapolation (linear in 
m^i) using smallest five masses. 



smallest quark masses is slightly higher than that ob- 
tained from the Xnn operator. For iViV-type field the 
mass difference is ~ 56(1) MeV and ~ 59(2) MeV for 
diquark-type interpolator. The positive mass splitting in 
both these channels could be interpreted as the evidence 
of repulsion. The chiral uncertainties are again estimated 
to be less than a percent. Similarly, positive mass dif- 
ferences that approximately remained constant with m~ 
and consistent with those obtained at our regular lat- 
tices was observed on a larger lattice volume of 6.24 3 fm 3 . 
Again the signature of repulsion at quark masses near the 
physical regime would imply no evidence of the resonance 
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FIG. 4: As in Fig. El but at a s = 0.26 fm. 

in the J PC = h channel. 

Finally, we performed a continuum extrapolation for 
the chirally extrapolated quantities in Fig. [5] Expecting 
that dominant part of scaling violation errors is largely 
eliminated by tadpole improvement, we adopt an a 2 s - lin- 
ear extrapolation for the continuum limit. We also per- 
form an a s -linear extrapolation to estimate systematic 
errors. Performing such extrapolations, we adopt the 
choice which shows the smoothest scaling behavior for 
the final values, and use others to estimate the system- 
atic errors. As can be seen from Fig. [5l the mass ratio 
again shows a weak dependence on the lattice spacing 
and varies only slightly over the fitting range. The four 
non-zero lattice spacing values of the ratio are within 0.01 
- 0.02 standard deviations of the extrapolated zero lattice 
spacing result. This will make for unambiguous and ac- 
curate continuum extrapolation. Using the physical nu- 
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FIG. 5: Continuum extrapolation of the mass ratios ex- 
tracted from diquark-type (solid circles) and NN (solid tri- 
angles) fields. The solid and the dashed curves are linear 
extrapolation to the continuum limit. 



cleon mass m p — 938 MeV, we obtain a continuum mass 
estimate of 1935(19) (60) and 1940(26) (72) MeV for the 
ground state hexa quark state from Xnn an d Xdq interpo- 
lators, respectively. The results from two extrapolations 
seem to be in good agreement, within errors. Given the 
fact that the ratio does not show any scaling violations, 
we could also quote the value of this quantity on our 
finest lattice, which has the smallest error. Nevertheless, 
order 4% errors on the finally quoted values are mostly 
due to the chiral and the continuum extrapolations. The 
quenching errors might be the largest source of uncer- 
tainty. Note however, that in the case of stable hadrons, 
this is not expected to be very important. It has been 
shown [l8j that with an appropriate definition of scale, 
the mass ratios of stable hadrons are described correctly 
by the quenched approximation on the 1 — 2% level. To 
this end we also calculated the pseudoscalar to vector 
meson ratio R$ p and pseudoscalar to nucleon mass ratio 
Rsn and found that in the continuum these ratios differ 
about 3% from their corresponding experimental values. 
So we quote our quenching errors to be less than five per- 
cent. This implies that dynamical quarks might play a 
less important role on the spectrum in question. 

Since the mass difference between the reported experi- 
mental pp mass and the physical 2m p continuum is ~ 30 
MeV, naively one may be tempted to interpret the results 
in Figs. and [Has a signature of the proton-antiproton 
bound state on the lattice. But the positive mass dif- 
ference observed in the range of pion mass differs from 
what is seen for the pp state studied on the lattice. The 
continuum results imply that the mass difference does 
not move into the continuum with an attractive interac- 
tion between p and p. This suggests that the observed 
signal is too heavy to be identified with the empirical 
pp(1859) and unlikely to be a signature of a possible 
proton-antiproton system as a six quark state. However 
the interaction between a baryon and an antibaryon can 
be as strong as that of baryon-baryon systems. Therefore 
we do not discard the interesting possibility that it could 
be a mixed molecular state of baryon and antibaryon 
bound by some long-range force of strong or electromag- 
netic nature or both. 

We have presented the results of the first lattice inves- 
tigation on the pp state employing improved gauge and 
fermion anisotropic actions, relatively light quark masses 
as well as smearing techniques to enhance the overlap 
with the ground state of the particle. Our analysis takes 
into account all possible uncertainties, such as statistical, 
finite-size, and quenching errors when performing the chi- 
ral and continuum extrapolations. On the basis of our 
lattice calculation we speculate that the state is not to 
be identified as a bound state of six quarks. However a 
thorough examination of this question would require the 
implementation of flavour SU(3) violation. The I = 0, pp 
state couples to Airr) [l(| through the ss component of the 
rj in the quenched approximation. By giving the strange 
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quark a larger mass would alter threshold which in turn 
would affect the manifestation of the bound state. As 
differences between meson masses in full and quenched 
QCD of order 60 MeV and 100 MeV or more for 
baryons have been observed [20] , one might wonder what 
effect the dynamics of full QCD could have on this state. 

Whereas the quenched approximation is found to give 
quite useful results at zero temperature, the fermions 
cannot be neglected if precise contact with nature is 
needed. One cannot yet rule out the possible existence 
of a pp bound state in full QCD. 

I would like to thank R. Jaffe, C. Michael, C. De- 
tar, S. Sasaki, C. Hamer and D. Leinweber for valu- 
able suggestions. It is a pleasure to acknowledge Shan 
Jin for his hospitality at CCAST and stimulating dis- 
cussions. We are grateful for the access to the com- 
puting facility at the Shenzhen University on 128 nodes 
of Deepsuper-21C. The computations were done using a 
modified version of the publicly available MILC code (see 
www.physics.indiana.edu/sg/milc.html). This work was 
supported by the Guangdong Provincial Ministry of Ed- 
ucation. 



[1] J.Z. Bai et al, Phys. Rev. Lett. 91 022001 (2003) 
[2] K. Abe et al, Phys. Rev. Lett. 88 181803 (2002) 



[3] K. Abe et al, Phys. Rev. Lett. 89 151802 (2002) 

[4] Particle Data Group, Phys. Rev. D 66, 010001 (2002) 

[5] E. Klempt, F. Bradamante, A. Martin and J.-M. Richard, 

Phys. Rep. 368, 119 (2002) 
[6] J.-M. Richard, Nucl. Phys. B (Proc. Suppl.) 86, 361 

(2000) 

[7] A. Datta and O'Donnell, Phys. Lett. B 567, 273 (2003) 

[8] J.L. Rosner, Phys. Rev. D 68, 014004 (2003) 

[9] Mu-Lin Yan, Si Li, Bin Wu and Bo-Qiang Ma, Phys. Rev. 

D 72, 034027 (2005) 
[10] Chong-Shou Gao and Shi-Lin Zhu, Commun.Theor.Phys. 

42, 844 (2004) 

[11] R.L. Jaffe, Phys. Rep. 407, 1 (2005); Nucl. Phys. (Prog. 

Suppl) 142, 343 (2005) 
[12] S. Sasaki, Phys. Rev. Lett. 93, 152001 (2004) 
[13] B. Lasscock, et al, Sasaki, Phys. Rev. D 72, 014502 

(2005) 

[14] C.J. Morningstar and M.J. Peardon, Phys. Rev. D 60, 
034509 (1999) 

[15] M. Okamoto, et al, [CP-PACS Collaboration], Phys. 

Rev. D 65, 094508 (2002) 
[16] P. Lacock, et al, [UKQCD Collaboration], Phys. Rev. D 

51, 6403 (1995) 
[17] D. Leinweber, A.W. Thomas, and R.D. Young, Phys. 

Rev. Lett. 92, 242002 (2004) 
[18] C. Gittringer et al., [BGR Collaboration] Nucl. Phys. B 

677, 3 (2004) 

[19] A. Khan et al., [PCAC Collaboration] Phys. Rev. D 65, 
054505 (2002) 

[20] R. Young, D. Leinweber, A. Thomas, and S. Wright, 
Phys. Rev. D 66, 094507 (2002) 



